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A density-matrix formalism which includes the effects of three-body ground- state correlations is 
applied to the standard Lipkin model. The reason to consider the complicated three-body correla- 
tions is that the truncation scheme of reduced density matrices up to the two-body level does not 
give satisfactory results to the standard Lipkin model. It is shown that inclusion of the three-body 
correlations drastically improves the properties of the ground states and excited states. It is pointed 
out that lack of mean-field effects in the standard Lipkin model enhances the relative importance 
of the three-body ground-state correlations. Formal aspects of the density-matrix formalism such 
as a relation to the variational principle and the stability condition of the ground state are also 
discussed. It is pointed out that the three-body ground-state correlations are necessary to satisfy 



o : 

(N . 

j , the stability condition. 

^ ■ PACS numbers: 21.60.Jz 



I. INTRODUCTION 

I ' Mean- field theories such as the Hartree-Fock theory (HF), the Haree-Fock Bogoliubov theory, the random-phase 
^ . approximation (RPA), and the quasi-particle RPA have extensively been used to study the ground states and collective 
^ ' excitations of atomic nuclei [1| . For more realistic theoretical treatment of nuclei such as inclusion of the ground-state 
P5 , correlations other than pairing correlations and the damping effects of collective excitations, however, we must go 
beyond the mean-field theories. The time-dependent density-matrix theory (TDDM) [2-4] which has been formulated 
I I by truncating the chain of the equations of motion for reduced density matrices up to the two-body level is one 
J> ' of such extended mean-field theories. It has been pointed out that a stationary solution of the TDDM equations 
\Q ' gives a correlated ground state and that the small amplitude limit of TDDM based on the correlated ground state 
, corresponds to an extended RPA (ERPA) including two-body transition amplitudes [5|. To test the reliability of the 

■ density-matrix formalism, we have applied it to solvable models [^-Q and found that the obtained results improve 
, on deficiencies of the mean-field theories ■ In the case of the standard Lipkin model @ where the interaction 

term contains only two particle - two hole excitations, however, we found that the ground states in TDDM become 
J slightly overbound as compared with the exact solutions. That the approximate ground states have lower energy than 

■ the exact ones contradicts the variational principle [l[ for the total wavefunction and should possibly be avoided, 
I ' although TDDM is not based on a Raleigh-Ritz variational principle. The aim of this paper is to clarify the origin 

^ . of such an unsatisfactory feature of TDDM in the standard Lipkin model. We will show that inclusion of three-body 
ground-state correlations, though it is complicated, drastically improves the results for the standard Lipkin model. 
^ , The paper is organized as follows: The formulation of TDDM with the three-body ground-state correlations jTll . [T2j 
and ERPA built on the TDDM ground state are given in sect. 2. Some formal aspects of TDDM and ERPA, a relation 
of the TDDM equations to the variational principle, the stability condition of the ground-state, and a comparison of 
ERPA with the self-consistent RPA (SCRPA) [Ij, [ij], which have not been pointed out in our earlier publications, 
are also discussed in sect. 2. The results for the standard Lipkin model are presented in sect. 3. Section 4 is devoted 
to the summary. 



II. FORMULATION 



Although the numerical calculations are performed for the Lipkin model, the formulation is presented using the 
following general hamiltonian 

H = J2mX')aiax' + l ^ {X^XMK>^2)a{A.''x',<^K^ (1) 

AA AiA2A^A2 

where t is the kinetic energy operator, w is a two-body interaction and a\^{ax) the creation (annihilation) operator of 
a nucleon in a single-particle state A. 
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A. TDDM 



In TDDM the ground state |0) is defined by the occupation matrix Uaa' ^ the two-body correlation matrix Ca^a'P' 
and the three-body correlation matrix Cap-^a' I3'i' given by 

naa' = {Q\ai,ac,\Q), (2) 
Caiia'P' = {Q\a^^,ayai3aa\0) - A{naa'ni3p,), (3) 

Capia'i3'-i' = (0|aJj,a]3,a^,a^a/3aQ|0) - y4,(naQ/ri/30/n^^/ -I- 5(nQ,Q'C/3^/3'^'))7 (4) 

where A and S mean that the products in the parentheses are properly antisymmetrized and symmetrized under 
the exchange of single-particle indices [3]. Equations of motion for Uaa' ^ Capa'P' and Cap-ya'P'j' can be obtained by 
truncating a coupled chain of equations of motion for reduced density matrices (the so-called Bogoliubov-Born-Green- 
Kirkwood-Yvon hierarchy) up to the three-body level jlj and can be written as 

ihj^n^^' = {0\[ala^,H]\0) = Fi{aa'), (5) 

ih^^Ca^pa^'P' = : {0\[alal,apaa,,H]\0) := F2ial3a'P'), (6) 

ih-^Cap^a'P'Y = : {0\[al^,a^f^,al^,ajapaa,H]\0) F3{al3-/a'l3'^'), (7) 

where : : means that the time derivatives of the second terms on the right-hand sides of eqs. ^ and (|3]) are 
subtracted. Since a four-body correlation matrix is neglected, the expectation values of four-body operators in eq. ([7]) 
are approximated by the products of riaa', Capa'P' and Cap-ya'P'^y'- The expressions for Fi, F2 and F^ are given in 
the Appendix, where the single-particle states which satisfy the HF-like mean-field equation 

/i(p)0,(l) = e„0„(l) (8) 

are used. Here, p is the one-body density matrix given by 1') = X]qq' "f^aa' 4'ai^)4'a' i^') ^^'^ numbers indicate 
spatial, spin and isospin coordinates. 

To obtain the ground state implies that all quantities Uaa' , Capa'P', Cap-ya'P'-y' and (pa are determined under the 
stationary conditions 

i^-^riaa' = Fi(aa') = 0, (9) 

irijCapc'P' = F2(a/3a'/3') = 0, (10) 

irij^Capya'P'-,' = FaiaMa'P'i) - 0. (11) 



and eq. ([5]). This task can be achieved using the gradient method [15|. The functional derivatives of eqs. (P))- ([TT|) are 
written as 



An \ / AFi 

AC2 = AF2 I = - I F2 I . (12) 




Sn SC2 
SF2 SF2 SF2 

Sn SC2 SC's . , , , 

SF3 Ma. Mk \ AC-i / \ AF- 

\ Sn SC2 SC3 

The inversion of the above matrix gives the following equation to be used in the gradient method 

n{N + l) \ ( n{N) \ / a c \ ( Fi{N) 



C2{N + l) = C2(N) \ -a\bde\ F2{N) | , (13) 
C^iN + l)) \C^{N) j {fghj \Fs{N) 

where n{N), C2{N) and C3{N) imply Uaa', Capa'P' and Cap-ya'P'-y' at the A^th iteration step, respectively, and 
a = SFi/Sn, b = SF2/Sn, c = SF1/SC2, d = SF2/SC2, e = SF2/SC3, f = SF^/Sn, g = SF3/SC2, and h = SF3/SC3. 
The expressions for these matrices are given in ref. |12| . The matrices a, b, d, /, g, and h depend on the iteration step 
N through n{N), C2{N) and C3{N). To solve eq. we start from a simple ground state such as the HF ground 

state where naa'i Capa'P' and Capya' p'y can be easily evaluated and iterate eq. (jl3p until convergence is achieved. 
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A small parameter a is introduced to regulate the convergence process. When the mean-field potential is present, eq. 

couples to eq. ^ through n^Q,/ . As was discussed in ref. [I^l , the matrix consisting of the functional derivatives 
of Fi , F2 and F3 on the right-hand side of eq. (fT3|) can be derived as the small amplitude limit of TDDM and has 
a close relation with the hamiltonian matrix of ERPA given in the next subsection [ll|, [l^ . This indicates that the 
ground state is not independent of the excited states in ERPA. Let us discuss this point in some more detail using 
the eigenstates of the following matrix equation, which corresponds to the small amplitude limit of eqs. ©-([Tl) [l2| . 




(14) 



Since the closure relation is given as [16 




E {£;y;s;)^i, (15) 



where (i* z*) is the left-hand eigenvector of eq. (ITil) and / is unit matrix, eq. (IT^ is written as 
n{N + l)\ f n{N) \ /x,\ fFi{N)\ 

C2{N+1) = CiN) \-aJ2—{ y, ] {x*^ y*, S;) F,iN) . (16) 



In eq. only the eigenstates with V,^ ^0 can contribute as is understood by multiplying eq. with (i* y* z*) . 
The above equation indicates that the occupation matrix and the correlation matrices are constructed from the 
eigenvectors of eq. (|14|) at each iteration step which have the same quantum numbers as the ground state: In the 
standard Lipkin model these states are two-phonon states. 

For the numerical calculations shown below we use not eq. but eq. because a considerable dimension 

size of the three-body part of eq. p4l) makes it impracticable to adopt eq. P^ which requires the eigenvalues and 
eigenvectors at each iteration step. In the matrix inversion of eq. (jl3p . however, mixing of unphysical states with 
ilp « is unavoidable. This problem can be removed by slightly shifting the unperturbed energies in o, d, and h by 
Ae so that the inverse can be taken. In the case of the Lipkin model the typical value of Ae used is 10~^e where e is 
the level spacing of the Lipkin model hamiltonian. The obtained results are independent of Ae because the product 

/F^iN)\ 

{i; y; s;) f,{n) (i?) 

V F^iN) J 

vanishes for the states with 17^ = 0. Since the inversion of the matrix on the right-hand side of eq. p^ is time 
consuming, the gradient method is supplemented with a time-dependent approach as will be explained below. We 
make a comparison with a simplified version of TDDM where the three-body correlation matrix Ca^^a' p'-y' is neglected. 
We refer this TDDM as to TDDMl. 



B. Extended RPA 

In this section we present our extended version of RPA (ERPA) [l^. We also discuss its close connection to TDDM. 
ERPA is formulated for the following excitation operator consisting of one-body and two-body operators 

= E ^AA' ■■ «I«A' : + E ^^lA.AiA^ 4i«LaA^«A; (18) 

AA AiA2A^A2 

where : : implies that uncorrelated parts consisting of lower-level operators are to be subtracted; for example, 

: a^j/flQ : = 0^^,00 , (19) 
■ al-ajj-a^aa : = a^a'^^p'^P'^a - AS{naa' ■ 0)3/0^ :) 

- [A{naa'npi3') + CaPa'P']- (20) 
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It is assumed that the operator Qj^ satisfies 



Q,.|*o> =0. 



(21) 
(22) 



Here, |^o) is the ground state in ERPA and |^^) is an excited state. The ERPA equations are obtained from the 
equations-of-motion method [l3l | 



(^'o|[[: al,al,apao. :, Qj,]|*o) = ^m(*o|[: al,al,af3ac, :,QJ,]|*o), 



(23) 
(24) 



where w^, is the excitation energy of 1^^^). In the evaluation of the matrix elements we approximate the ERPA ground 
state l^o) by |0) which satisfies eqs. ([5|)- pT|) . The ERPA equations can then be written in matrix form 



^1 Ti 

T2 5*2 



where the matrix elements are given by 



A{aa' : XX') 


= (0 


B{al3a'l3' : XX') 


= (0 


C{aa' : XiX2X'iX'2) 


- (0 


D{aPa'P' : AiAaA'^A^) 


= (0 


Si{aa' : AA') 


- (0 


Ti(aa' : AiAaA'^A^) 


- (0 


T2{aPa'f3' : AA') 


= (0 


S2{a(3a'P' : XM'^X'^) 


= (0 



[: a^.tta ■■,H],: a{ax> :]|0), 

[: al^.alj.apaa ■.,H],: alax> :]|0), 

[: al,aa -^H],: a{^a{^ax'^ax'^ :]\0), 

[: al^,a''p,ai3aa ■,H],: a{y^^ax'^ax'^ :] 

: al,aa :, : a^OA' :]|0), 

: al,aa : a{^a{^ax'^ax'^ :]|0), 

■ ai-apa^fla : a\ax' :]|0), 

: a^,ajj,a^aa :, : al^al^ax'^ax'^ :]|0). 



(25) 



(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



At this point we want to insist on the fact that in eqs.(l^- pT|) antisymmetrisation is fully respected and, thus, the 
matrix elements entering eq. (1251) also respect antisymmetrisation fully. This is at variance of most other extensions 
of the RPA approach. 

In the following we show that eq. (|25l) can also be derived from the small amplitude limit of TDDM (eq. (ITil) V We 
transform the eigenvector in eq. (|14p using an extended norm matrix including three-body components as 




where 



Si Ti T13 
T2 S2 T23 
T31 T32 T33 



T31 = (0|[: 0^^,0)3, ajj,,a^a;3aa :,: a^flA' :]|0), 

732 = {0\[: al,al,al^,aja(3aa a{y^^ax'^ax'^ ■.]\0) 

Ti3 = {0\[: al,aa a{^a{^_al^ax'^ax'^ax'^ --M, 

T23 = {0\[: al,al,af!aa al^al^al^ay^ax'^ay^ --M, 

T33 = {0\[: al,al,al,a^ai3aa :,: a{^a{^a{^ax'^ax'^ax'^ 



Then eq. (|14l) is written as 

aSi + CT2 



aTi + CS2 aTi3 + CT23 

bSi + dT2 + eT^i bTi + 6,82 + eT'32 hTi^ + dT23 + er33 
fSi + gT2 + hTzi fTi + 56*2 + hT32 /T'13 + .gT'23 + /1T33 




Si Ti Ti3 
T2 S2 T23 
T31 T32 T33 




(34) 



(35) 
(36) 
(37) 
(38) 
(39) 



(40) 



The matrix elements on the left-hand side of the above equation can be written in the form of the double commutators 
with the hamiltonian [l^ except for those with /, g and h: Since [: a^^iO^i^ial^ia^apaa '■, H] contains four-body operators. 
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the matrices fSi + 5T2 + hT^i, fTi + gS2 + /1T32, and /T13 + (7T23 + hT^s cannot be expressed using the double 
commutator with the hamiltonian. The matrices A, B, C, and D in eq. (j25p have the fohowing relations with a, b, 
c, d, and e A = aSi + 0X2, B = bSi + dT2 + eTji, C = aTi + cS'2, and 13 = bTi + (iS'2 + eTi2- The expressions 
for iSi, Ti, T2, 5*2, Tsi, and T32 are given in ref. 12]. Thus it is shown that eq. ([^5)) is obtained from eq. (|40p by 
neglecting the three-body sections. _ 

As has been discussed in detail in ref. the hamiltonian matrix on the left hand side of eq. (|25p is hermitian 
due to the ground-state conditions eqs. ©-([H]) which guarantee the Jacobi's identity 

where P and Q are either one-body or two-body operators. It is interesting to note that reversing the order of 
arguments and imposing hermiticity of the matrix on the left-hand side of eg. ipSj) one can arrive at eqs. (P|)-([TT|). We 
remind in this context that Rowe |13| achieved hermiticity of his equation of motion method taking the arithmetic 
mean of the off diagonal elements. This somewhat artificial procedure finds a natural solution with the considerations 
given in our present procedure. It, therefore, can be stated that in the equation of motion method the non hermitian 
part of the matrices has to be put to zero what imposes the fulfillment of extra equations. 

As mentioned above, the hamiltonian matrix of eq. (1401) is not hermitian because of the asymmetry in the three- 
body parts: It can be stated that eq. (j25p is obtained omitting the non- hermitian parts of eq. (j40p . The commutator 
[P, Q] for three-body operators P and Q contains three-body, four-body and five-body operators. In order to fulfill 
the condition eq. (PT|) and make an extended RPA with the three-body amplitude hermitian, therefore, we need 
to consider two additional ground-state conditions for four-body and five-body operators. The problem of a further 
extended RPA with the three-body amplitude is beyond the scope of this paper. 

The ortho- normal condition of ERPA (eq. ((25]) ) is given as [19| 

Xn{% |)(i:')-±V, (42) 

where the negative sign is for a negative-energy solution. Equation ([^5]) is equivalent to the second RPA [13, ITsj 
when the ground-state correlations are neglected: Neglect of the ground-state correlations means that riaa' = 5aa' (0) 
for occupied (unoccupied) states, Capa'p' ~ and Cap-ya'P'-y = 0. In this limit, the one-body section of eq. (|25p . 
Ax^ = ujfj.Six^, is equivalent to standard RPA. We make a comparison with a simplified version of ERPA where the 
three-body correlation matrix Cap-ya' p'-y' is neglected in the calculation of the ground state. We refer this ERPA as to 
ERPAl. When eq. ((25|) is solved numerically, first the norm matrix on the right-hand side of eq. ((25|) is diagonalized. 
Then eq. (|25p is solved in the truncated space consisting of the eigenstates of the norm matrix with nonvanishing 
eigenvalues. 



C. Stability condition 

We discuss the relation of the ground state conditions eqs. ([9]) and (jlOp to the variational principle for the energy. 
We also discuss the stability condition of the ground state. Suppose the following variational energy 

E= {Q\e-'^He'^\Q), (43) 

where F is an arbitrary hermitian operator given as 

F = ^ faa' ■ a-iaa' : -f- ^ Fafja'p' ■ ala^ja/H'fla' : ■ (44) 

act' ol^ol' (5' 

The energy E can be expressed as 

E = {0\H\0)+i{0\[H,F]\0} + ^{0\[[F,H],F]\0} + ---- 
= {0\H\0)+iY,{0\[H,:aiao.' :]|0)/„„. 

aa' 

+ i ^ {0\[H,: al^a''^ai3>aa' ■■]\0)Fapa'i3' 

a0a'P' 



+ 2^^* ^*nB D [f)+----- (45) 
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The stationary conditions eqs. ([9]) and (ITOl) correspond to 5E/5faa' = and 5E/5Fapa'P' = 0, respectively, and 
the single-particle hamiltonian of eq. ([8|) is obtained from the variation 5{Q\H\Q) / 5naa' = haa' ■ The condition eq. 
([TT|) can also be obtained by adding a three-body operator to eq. (03]). If the total wavefunction |0) were used to 
calculate riaa', Capa'ii', Caj3-ya' p'-y' , and also the four-body correlation matrix, eas.([9])- (fTT|) would correspond to the 
variational principle for the total wavefuntion. However, it is impracticable to find |0) itself. In TDDM, eqs.(|9|)-(fTT|) 
are used to obtain not the total wavefunction but Uaa' , Gapa' p' , and Cap-ya' p'-y' by approximating a four-body density 
matrix with lower level density matrices. In this sense, TDDM deviates from the variational principle for the total 
wavefunction, although eqs.(|9|)- (llll) can be formally derived from the variational principle. 

The stability matrix in the last line on the right-hand side of eq. (j45p is nothing but the hamiltonian matrix of 
eq. (PSI) . As was mentioned above, the condition for the three-body matrix eq. ([TT|) in addition to eqs. © and ([TU]) 
is necessary to make the stability matrix hcrmitian although the operators eqs. ([T8| and (|44p are at most two-body 
ones. The stability condition that |0) corresponds to a minimum in the energy surface is given by 

(/* ^)[i)>0 (46) 

for arbitrary faa' and Fapa'P'- This condition is satisfied when the stability matrix is positive definite. In this case, 
eq. (|25p has only real eigenvalues since the norm matrix in eq. (|25p is hermitian |Jj]. 

D. Relation to self-consistent RPA 

In this subsection we discuss some relation of the one-body part of eq. ([25)1 

Axf" ^ ujf.Six" (47) 
to SCRPA [14]. The matrix A is explicitly written as 

A{aa' : XX') = (Ca - ^a'){n\>a'5a\ ~ na\5a'\') 

77' 

- ({-f'^\v\a'X)An\'^ - (7'A'|u|a'7)A%A)"-Q7'] 

- X! (("7|^'|7'7")C'7'7"A7<5q'A' + (77'I«I"'7")C'a'7"77'^qa) 
77'7" 

+ {X'-f\v\a''y')A 

- ^((aA>|77')C^7'a'A + (77>|a'A)C„A'77'), (48) 
77' 

where the subscript A means that the corresponding matrix is antisymmetrized. The norm matrix 5*1 is given by 

Si{aa' : AA') nya'SaX - "-oA^a'A'- (49) 

The first two terms with Capa' p' on the right-hand side of eq. (|48)) describe the self-energy of the particle - hole state 
due to ground-state correlations pol |. The last four terms with Capa'P' may be interpreted as the modification of the 
particle - hole interaction caused by ground-state correlations [lO] • Equation (|47ll is formally the same as the SCRPA 
equation. Both equations include the effects of ground-state correlations through Uaa' and Capa' p' ■ The difference 
between eq. (1471) and the SCRPA equation lies in the fact that Uaa' and Capa' p' are self-consistently generated from 
^aa' SCRPA, while they are given by eqs. (O and ((TU| , independently of eq. (|T7l) . Another difference stems from 
the fact that eq. © is used in SCRPA to determine the optimal single particle basis whereas the occupation numbers 
are obtained from a separate relation [l3 |. 

III. APPLICATION TO LIPKIN MODEL 
A. Lipkin model 

The Lipkin model ^ describes an N-fermions system with two N-fold degenerate levels with energies e/2 and — e/2, 
respectively. The upper and lower levels are labeled by quantum number p and —p, respectively, with p = 1, 2, N . 
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We consider the standard hamiltonian 

H^eJ, + ^{Jl + Jl), (50) 

where the operators are given as 

a_p^a_p), (51) 

N 

(52) 

p=l 

Since eq. dH) keeps the property pp — 0, the single-particle state in eq. ([8|) is equivalent to the original 

single-particle basis labeled by —p and p. When the matrix inversion in eq. (|13p is taken, it is essential to preserve 
symmetry properties of all the matrices. Therefore, we use the so-called m scheme to define each matrix throughout 
our numerical calculations: For example, n_j,_p, n^pi^pi, C^p-p'pp', C^p'^ppp', C-^p^p'p'p and C^p'^pp'p with p and 
p' = 1,2, ...,A^ are treated as independent quantities. Of course, eqs. © and ([TU| guarantee n_p_p=n_p'_p' and 
C-p-p'ppi — —C-pi-pppi = —C-p-p'p'p — C-pi-ppip. Let us mention that due to the use of the m-scheme the 
dimension of the matrices can become considerable, even for the case of the Lipkin model. One is forced to do that, 
otherwise the full antisymmetry can not be maintained. 



1 ^ 



J+ = Jl 



B. Ground state 



Since eq. (I13p involves the time-consuming inversion of a large matrix, we first find an approximate solution for 
eqs. @ - ([TTt to be used in the gradient method using the following time-dependent approach [H, [2l|: Starting 
with the HF ground state, we solve eqs. dS])-© using the time-dependent interaction V{t) — V x t/r with large r 
[21I [2^ . We tested the reliability of this time-dependent method for an iV = 2 system, to which the truncation of 
the reduced density matrices up to the two-body level should give the exact solution. Since there are no three-body 
density matrices, Eq. ^ for the N — 2 system is modified to 

d 

+ ^ [(q;/3|i;|AiA2)paiA2q'/J' - (AiA2|f;|a';3')Pa/3AiA2], (53) 

Ai A2 

where PajSa'P' = ■Ainaa'npfji) -\- Cafja'iS'- We found that r = 27r?i/e x 10 is sufficiently large to suppress the spurious 
mixing of excited states and the obtained ground-state energy is equivalent to the exact one — Ve^ + within 
numerical accuracy. After the time-dependent calculation eq. (|13|) is then solved for a few hundred iterations to 
achieve the conditions eqs. © - (jlip for each matrix element, which guarantee the hermiticity of the hamiltonian 
matrix of eq. (1^ . Since the time dependent approach already gives a good approximate solution for eqs. (O - pD) . 
the change in the total energy during the iteration process of eq. (jTS)) is negligible. 

The ground-state energy in TDDM (solid line) is shown in fig. [l]as a function of x = [N — ^)V/e for = 4. The 
results with TDDMl and the exact ground-state energies are shown in fig. [1] with the dotted and dot-dashed lines, 
respectively. As seen in fig. [U the ground-state energies in TDDMl are overestimated. This unpleasant feature of 
TDDMl is removed by the inclusion of the three-body correlation matrix and the agreement with the exact solution 
is much improved. In order to see the ground-state properties in more detail, we show the occupation probability 
of the upper level npp and the two-body correlation matrix Cppi-p-pi in figs. [5]and|31 respectively. The solid and 
dotted lines denote the results in TDDM and TDDMl, respectively. The exact values are shown with the dot-dashed 
lines. Figures [5] and [3] show that the agreement of these matrices with the exact ones is also improved by the inclusion 
of the three-body correlation matrices. As is understood from eqs. \A2\ and (|A7I) in the Appendix, the three-body 
correlation matrix interferes in particle - particle and particle - hole correlations and presumably plays a role in 
screening these two-body correlations. In order to investigate the effects of the three-body correlations in larger N 
systems, we calculate the ground-state energies in TDDM for N = 8. As shown in fig. 21 an improvement similar to 
the = 4 case is achieved. 

The importance of the three-body correlation matrix in the standard Lipkin model is strongly in contrast to the 
case of an extended Lipkin- model hamiltonian which we have previously studied [lOl , [l^ . The extended Lipkin model 
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FIG. 1: Ground-state energy in TDDM (solid line) as a function of x = (A^ - l)|V'|/e for A*' = 4. The dotted line depicts the 
results in TDDMl where the three-body correlation matrix is neglected, and the dot-dashed line the exact values. 
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0.2- 




FIG. 2: Occupation probability of the upper state as a function of x for A'^ = 4. The meaning of the three lines is the same as 
in fig. [U 
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FIG. 3: Two-body correlation matrix Cpp'_p_p/ as a function of x for A'^ = 4. The meaning of the three lines is the same as in 

fig.m 
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FIG. 4: Ground-state energy as a function of x for N — 8. The meaning of the three hues is the same as in fig. [T] 
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FIG. 5: Strength function of the one-phonon state in ERPA (sohd line) for x — 1 A*' = 4. The dotted and dot-dashed hues 
depict the resuh in ERPAl and the exact solution, respectively. 
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FIG. 6: Strength function of the one-phonon state for x ~ 1-5 and = 4. The meaning of the three lines is the same as in fig. 
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FIG. 7: Excitation energies of the one-phonon state (lower part) and the two-phonon state (upper part) as a function of x for 
= 4. The filled and open squares depict the results in ERPA and ERPAl, respectively. The filled and open triangles indicate 
the results in IRPA and IRPAl, respectively. The dashed line depicts the results in RPA. The exact solutions are given by the 
solid lines. 



has the additional particle scattering term {U term) 0, Q 

if/[J,(J+ + J_) + (J+ + J_)J,]. (54) 

This U term generates a mean-field potential [8] and gives nonzero ofF-diagonal elements of the occupation matrix, 
that is, rip_p ^ and n„pp ^ 0. Therefore, a significant portion of the interaction energy including the V term in eq. 
(1501) is carried as the mean-field energy jlO|, almost independently of the strength of the U term. As a consequence, 
the two-body correlation matrix and the correlation energy become quite small in the extended Lipkin model. The 
effect of the three-body correlations on the ground-state energy was found even smaller when the U term was included 
[l^ . In the case of the extended Lipkin model the ground-state energies calculated in TDDMl are always above the 
exact values and such an overbound problem as in the standard Lipkin model does not occur. Now the importance of 
the three-body ground-state correlations in the standard Lipkin model is understood as a consequence of lack of the 
mean-field energy, which makes the truncation of the reduced density matrices up to the two-body level unreliable 
and enhances the relative importance of the three-body correlations. 

The standard Lipkin model has deformed HF solutions with n_pp ^ for x > 1 which consequently carry 
the mean-field energy which sums up already a lot of correlations. One may then think that the truncation up to 
the two-body level be reliable when the two-body correlation matrix is defined using the deformed HF single-particle 
states. To answer this question, we performed a simplified gradient-method calculation where the three-body parts in 
eq. ([T^ are neglected and a deformed HF state at x = 1-5 is used as the initial state. We found that the converged 
results are the same as those shown in figs. [T][3l The residual interaction plays a role in restoring the symmetry of eq. 
(1501) which is violated by the deformed HF solution. In other words starting with a 'deformed' solution, at the end of 
the iteration cycle the solution is back to 'sphericity'. Thus the use of the deformed HF basis does not improve the 
two-body level approximation for the standard Lipkin model. This somewhat surprising feature may be due to the 
small number of particles considered. It is known that standard deformed RPA becomes exact in the macroscopic 
limit. Therefore, we suppose that considering higher particle numbers, the deformation will catch on when we start 
with the deformed HF state. 



C. Excited states 



For the evaluation of the excited states, we use the ERPA equation (25). First we present the results for the 
one-phonon states in the N — A system excited by the operator Q = J+ + J^. Figures [5] and [S] show the strength 
functions of the one-phonon states at x = 1 and x = 1.5, respectively. To facilitate easy comparison among various 
calculations, we smooth the strength functions with an artificial width rpwHivi/e = 0.025. The solid and dotted lines 
denote the results in ERPA and ERPAl, respectively. The exact results are shown with the dot-dashed lines. The 
results in ERPA almost coincide with the exact ones, while those in ERPAl are located above the exact solutions. 
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FIG. 8: Strength function of the two-phonon state in ERPA (solid line) for x — 1 ^-nd A'' = 4. The dotted and dot-dashed lines 
depict the result in ERPAl and the exact solution, respectively. 
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FIG. 9: Strength function of the two-phonon state for x ~ 1-5 and A'^ = 4. The meaning of the three lines is the same as in fig. 

m 



Note that RPA breaks down at x = 1 m as shown in fig. H At x = 1.5 the effect of the three-body correlations on the 
one-phonon state is drastic because it significantly improves the occupation probabilities and the two-body correlation 
matrices as shown in figs [2] and [3] Thus it is found that the inclusion of the three-body ground-state correlations also 
gives a better description of the onc-phonon state in the standard Lipkin model. The x dependence of the excitation 
energy of the one-phonon state is shown in fig. [7] for various approximations. The filled triangles indicate the results 
obtained from eq. (|47p using the same Uaa' and Capa'iS' as those used in ERPA. This approximation is referred to 
as the improved RPA (IRPA). Similarly, the results obtained from eq. (|47)) with the same Uaa' and Capa'f)' as those 
used in ERPAl are referred to as the IRPAl results (open triangles). The excitation energies of the one-phonon state 
in IRPA and IRPAl are larger than the exact values and increase with increasing x- This is explained by the increase 
in the self-energy terms in eq. pS)) . Since the deviation of the occupation matrix and correlation matrix from the 
exact values is larger in TDDMl than in TDDM, the excitation energies in IRPAl are larger than those in IRPA. The 
difference between the results in IRPA and ERPA and also between IRPAl and ERPAl is due to the coupling of the 
one-body amplitudes to the two-body ones. Since the parity of the number of particle - hole pairs is conserved in the 
standard Lipkin model, the one-body amplitudes can couple to the two-body amplitudes that have the same parity: 
For example, Xp_p couples to X'^^, , which expresses excitations built on the two particle - two hole configurations 
in the ground state. As shown in fig. [71 the coupling to the two-body amplitudes is essential to obtain the one-phonon 
states with accurate excitation energies. 

Now we discuss the two-phonon states excited by = ( J+ + j_)^. The results at x = 1 and 1.5 are shown in 
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figs. [8] and [9l respectively. Due to the parity conservation, the two-phonon state excited by does not couple to 
the one-phonon state in the standard Lipkin model. The excitation energies both in ERPA and ERPAl are higher 
than the exact value and the improvement of the two-phonon state due to the inclusion of the three-body correlations 
is not so large as in the case of the one-phonon state. This indicates that the correlations among the two-body 
amplitudes are not sufficient to lower the energy of the two-phonon state, suggesting the importance of the coupling 
to the three-body amplitudes which are neglected in ERPA and ERPAl. The x dependence of the excitation energy 
of the two-phonon state is also shown in fig. [T] 

We also studied the stability condition by diagonalizing the stability matrix eq. (|46p . We found that all the 
eigenvalues associated with the physical operators such as J+, J_ and which consist of the hamiltonian eq. ((50)) 
are positive definite. This corresponds to the real eigenvalues of eq. (j25p for the one-phonon and two-phonon states 
excited by Q = J+ + J- and . However, the stability matrix has some negative eigenvalues for other degrees of 
freedom, though their absolute values are small; They are at most 1% of the maximum positive eigenvalue for the 
physical operators. This means that the stability condition is not always satisfied for unphysical states. These states 
correspond to non-collective phonons as discussed in ref. [131 ■ 



IV. SUMMARY 



The density-matrix formalism which includes the effects of three-body correlations on the ground state was applied 
to the standard Lipkin model. It was found that the inclusion of the three-body correlations removes unpleasant 
features of the two-body level approximation and drastically improves the ground-state properties. It was discussed 
that the importance of the three-body correlations in the standard Lipkin model is attributed to the fact that it 
does not have the mean-field contributions when the conservation of the number of particle - hole pairs is respected. 
The extended RPA built on the ground state with three-body correlations was applied to the one-phonon and two- 
phonon states. It was found that the spectrum of the one-phonon state is drastically improved by the inclusion of 
the three-body correlations. In the case of the two-phonon states, however, the agreement with the exact solutions is 
not so good as in the case of the one-phonon states, which suggests the importance of the coupling to the three-body 
amplitudes. Inclusion of the three-body correlations for realistic nuclei may be impracticable. Fortunately, mean- 
field effects are quite important in nuclei and the mean-field theories give good first description for ground states 
and collective excitations. Therefore, the truncation up to the two-body level may be justified as in the case of the 
extended Lipkin model. When the three-body correlation matrix is neglected, the hermiticity of the hamiltonian 
matrix in ERPA is not guaranteed. Our applications so far performed for realistic cases indicate that this does not 
cause any serious problems. 



Appendix A 

Fi, F2 and F3 in eqs. (in))- (|TT|) are shown. The single-particle states which satisfy the HF-like equation hi^a — ^4>a 
are used. 

Fi{aa') = (ea - ea')naa' + ^ iC\^\2a'xA<^^3\v\>^l>^2) ~ Ca\3XiX2{^1^2\v\a' X3)), (Al) 

Ai A2 A3 

F2 («/?«'/?') = (e^ + £0 — ea' — <il3')Cal3a'l3' + Bapa'l3' + Pafia' P' + Hafja'P' + Taf3a'P', (A2) 

F3{al3ja'l3'y) = (e^ + ejs + - e^' - ep- - e^,)Cai3^a' p'-i' 

- r{a'P'iap-i) + r{f3'a'-f'af3j) + r{j'P'a'al3j) 
+ J{al3-ia'l3'i) + J{l3-faa' P'-/') - J{a-/l3a' /S'-f') 

- J* {a'P'iap-i) - J* {f3'j'a'al3j) + J* {a'-f' /S'ap-f) , (A3) 

where 

BaPa'P' = ^ (AiA2|w|A3A4)a[((5qAi - ?^aAl)('5^^A2 - "/3A2)"-A3a'"-A4/3' 
Al A2 A3A4 

- na\^ni3x^{5x^a' - nx^a'){5\ii3' - nx^ii')], (A4) 

Papa'fi' — ^ (A1A2IWIA3A4) [(5qAj(5^A2 ^ ^qAi»^^JA2 — "-QAi^/3A2)C'A3A4a'/3' 
Al A2A3A4 
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— {Sxsa'Sx^p' — S\^a'n\^i3' — riAga' (5a40' Aa ] j (A5) 
Ai A2A3A4 

— <5/3'A4(«/3A2C'a3qAiq' — "qA2Ca3;9Aiq')]i (^6) 

Tapa'fi' = ^ (AiA2|w|A3A4)[(5QAiCA3A4/3a'A2;9' + '^/3A2CA4A3Q/3'Aia' 
Ai A2A3A4 

— (5Q'A3C'aA4^AiA2^' ~ '^^'A4C'/3A3qA2Aiq'], (A7) 

I{aPja' f3'j') ^ ^ {(aA3|i;|AiA2)(n7A3C'AiA2/3a'^'7' - 'T-0A3C'aiA27q'/3'7' + C'AiA2a'/3'C'^77'A3 
Ai A2A3 

~ C'Ai A2a'7'C'/37^'A3 ~^ CxiX20' j' C-fija' X3) 

+ (aA3|^'|AlA2)A[^^AlQ'C'A2070'7'A3 — nXil3'Cx2l3^a'^'X3 + «Ai7'C'a2/37q'/3'A3 

^~ C*A2/37'A3 C*Ai7a '/3' + C'a2 'A3C'a 1^7'^' + C'A27/3'A3CA4^a'7' Ca277'A3 CAi^ct'/3' 

^ C'a2/9q'A3Cai77'/3' ~ C'a2/3^'A3Cai7q'7' 

+ njX3inXia'Cx2l3p'j' — nXii3'Cx2l3a'',' + "-Ai7'C'a2/3q'^') 

~ 'T'^A3 ('T'Aiq'Ca27^'7' ~ ^Ai^'C'a27q'7' + 7^Ai7' Ca27q'/3' ) 

+ ?T-Aiq'(«A2/3'C'/377'A3 — "A27'C'/37^'A3) + «Ai^'«A27'C'/37q'A3]}, (A8) 

J(a/37a'/3'y) - ^ [(a/3|t;|Ai A2)A(nA4a'CA 27^'7' — "Ai/3'C'a2 CA27a'/3') 
Ai A2 

+ (a/?|i;|AiA2)CAiA27a'/3'7']- (A9) 

The matrix Bafia' does not contain Capa'/i' and may be called the Born term, whereas Papa' I3' and Ha^a'ti' contain 
Capa'p' and describe particle-particle (and hole-hole) and particle-hole correlations to infinite order, respectively. 
The last term on the right-hand side of eq. (|A2p expresses the contribution of the three-body correlation matrix. 
The matrix J{af3^a' fi'^') includes particle-particle (and hole-hole) correlations, while I{a/3'ja' /S'-f') contains both 
particle-particle and particle-hole correlations. 
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